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ABSTRACT 

We have deduced the cosmic microwave background (CMB) temperature in 
the Coma cluster (A1656, z = 0.0231), and in A2163 iz = 0.203) from spectral 
measurements of the Sunyaev-Zel'dovich (SZ) effect over four passbands at radio 
and microwave frequencies. The resulting temperatures at these redshifts are 
Tcoma = 2.789±g;g|g K and T A2W3 = 3.377^;^ K, respectively. These values 
confirm the expected relation T(z) = Tq(1 + z), where To = 2.725 ± 0.002 K 
is the value measured by the COBE/FIRAS experiment. Alternative scaling 
relations that are conjectured in non-standard cosmologies can be constrained by 
the data; for example, if T(z) = T (l + z) l ~ a or T(z) = T [l + (1 + d)z], then 
a = — 0.16loj$2 an d d = 0.17±0.36 (at 95% confidence). We briefly discuss future 
prospects for more precise SZ measurements of T(z) at higher redshifts. 



Subject headings: cosmic microwave background — cosmology:observations — galax- 
ies:clusters:individual(Al656,A2163) 
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1. Introduction 

The (present) CMB temperature was precisely measured by the Far Infrared Absolute 
Spectrometer (FIRAS) on board the COBE satellite, To = 2.725 ±0.002 K, in the frequency 
range 2-20 cm -1 (Mather et al. 2001). These measurements essentially rule out all cosmo- 
logical models in which the CMB spectrum is non-Planckian at z — 0. Models with a purely 
blackbody spectrum but with a different T(z) dependence than in the standard model are, 
however, unconstrained by the FIRAS database. Also unconstrained are models with spec- 
tral distortions that are now negligible, but were appreciable in the past. A specific example 
is the relation T(z) = T (l + z) l ~~ a , where a is a parameter of the theory (see, e.g., Lima et 
al. 2000). More generally, models in which ratios of some of the fundamental constants vary 
over cosmological time are also of considerable interest. 

So far, T(z) has been determined mainly from measurements of microwave transitions 
in interstellar clouds which contain atoms and molecules that are excited by the CMB (as 
reviewed by LoSecco et al. 2001). The temperature has been determined in the Galaxy, as 
well as in clouds at redshifts up to z ~ 3 (Levshakov et al. 2001). These measurements 
are affected by substantial systematic uncertainties stemming from the unknown physical 
conditions in the absorbing clouds (Combes and Wiklind 1999). 

The possibility of determining T(z) from measurements of the Sunyaev-Zel'dovich (SZ) 
effect had been suggested long ago (Fabbri, Melchiorri & Natale 1978, Rephaeli 1980). (For 
general reviews of the effect and its cosmological significance, see Rephaeli 1995a, Birkinshaw 
1999.) The proposed method is based on the steep frequency dependence of the change 
in the CMB spectral intensity, Al, due to the effect, and the weak dependence of ratios 
Al (vi) / Al (v>j) of intensity changes measured at two frequencies (i/j, Vj) on properties of the 
cluster (Rephaeli 1980). Because of this, and the fact that - in the standard cosmological 
model - the effect is essentially independent of z, SZ measurements have the potential of 
yielding much more precise values of T(z) than can be obtained from ratios of atomic and 
molecular lines. With the improved capability of reasonably precise spectral measurements 
of the SZ effect, the method can now be used to measure T(z) in nearby and moderately 
distant clusters. Here we report first results from spectral analysis of SZ measurements in 
the Coma and A2163 clusters of galaxies. 
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2. T(z) from SZ 

The CMB intensity change due to Compton scattering in a cluster can be written in 
the form 

' ' ! ' ' Oh{x)-P + R(x,0,P) 



A/ h 2 c 2 (e* - l) 2 J dT 

where x = hvjkT is the non-dimensional frequency, 9 = kT e /mc 2 is the electron tem- 
perature in units of the electron rest energy, and (5 is the line of sight (los) component 
of the cluster (peculiar) velocity in the CMB frame in units of c. The integral is over 
the Compton optical depth, r. Both the thermal (Sunyaev & Zeldovich 1972) and kine- 
matic (Sunyaev & Zel'dovich 1980) components of the effect are included in equation (1), 
separately in the first two (additive) terms, and jointly in the function R(x,9,/3). In the 
non-relativistic limit (which is valid only at low electron temperatures and frequencies) the 
spectral dependence of AI is fully contained in the product of the x-dependent pre-factor 
times the function fi(x) = x(e x + l)/(e x — 1) — 4. The more exact treatment of Compton 
scattering in clusters necessitates a relativistic calculation (Rephaeli 1995b) due to the high 
electron velocities. The function R(x,9,(3) includes the additional spectral, temperature, 
and (cluster) velocity dependence that is obtained in a relativistic treatment. This function 
can be approximated by an analytic expression that includes terms to orders 9 5 and 2 9: 



(2) 



R(x, 9, (3) ~ 9 2 [f 2 (x) +9f 3 (x) +9 2 f 4 (x) +9 3 f 5 (x)\ -(59 [ 9l (x)+9g 2 (x)\ +(3 2 [l + 9g 3 (x) 

The spectral functions /j and gi were determined by Itoh et al. (1998), Itoh et al. (2002), and 
Shimon & Rephaeli (2002). For our purposes here this analytic approximation is sufficiently 
exact even close to the crossover frequency. The non-relativistic limit, R = 0, applies if the 
sum of all these terms can be ignored at the desired level of accuracy. 

The z dependence of AI is fully determined by the functions v = v(z), and T = T(z). 
The temperature-redshift relation may assume various forms in non-standard cosmologies; 
here we consider two examples. In the first, T(z) = T(0)(1 + z)^~ a \ where a is taken to 
be a free parameter, but with the standard scaling v = z/ (l + z) unchanged. With these 
relations the non-dimensional frequency obviously depends on z, x = x (l + z) a , if a ^ 0; 
here, x = hu o /kT(0). Another functional form which seems also to be of some theoretical 
interest is T(z) = T(0)[l+(l+d)z] (LoSecco et al. 2001), for which x = x (l+z)/[l+(l+d)z}. 
Obviously, in the standard model a — d — 0. 

For a slow moving (f3 < 10~ 3 ) cluster, the expression for AI in the non-relativistic 
limit depends linearly on the Comptonization parameter, y = J 9dr, which includes all 
dependence on the cluster properties. A ratio of values of AI at two frequencies is then 
essentially independent of these cluster properties. In the more general case, the first term 
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in the square parentheses in eq. (1) still dominates over the other two, except near the 
crossover frequency (whose value generally depends on T e , except in the non-relativistic 
limit where x c = 3.83; Rephaeli 1995b, Nozawa et al. 1998, Shimon & Rephaeli 2002), where 
the sum of the temperature dependent terms vanish. For values of x outside some range 
(roughly, 3.5 < x < 4.5), the dependence of AI, and - particularly, a ratio of values of 
AI - on j3 is very weak since the observed temperature range in clusters corresponds to 
0.006 < 6 < 0.03, whereas typically f3 < 0.002. 



3. Data Analysis 

We have analyzed results of SZ measurements in the Coma cluster (A1656) and A2163. 
Measurements of Coma, z = 0.0231 ± 0.0017 (Struble and Rood 1999) were made with the 
MITO (DePetris et al. 2002) telescope in 20 hours of integration. We also use the result of 
measurements at 32 GHz made with the OVRO 5.5m telescope (Herbig et al. 1995, Mason 
et al. 2001). A2163, z = 0.203 ± 0.002 (Arnaud et al. 1992), was observed with the SuZIE 
array (Holzapfel et al. 1997a), and with both the OVRO and BIMA interferometric arrays 
(LaRoque et al. 2000). 

When observing a cluster the SZ part of the measured signal is 

POO 

ASi = GiAVt\i / A/(0) ei (i/)di/, (3) 
Jo 

where Gi is the responsivity of the i th photometric channel, AQ\i is the corresponding 
throughput, and €i(v) is the spectral efficiency The full measured signal includes also 
contributions from the atmosphere, CMB anisotropies, and - at very high frequencies - 
also emission from dust. Multifrequency observations allow us to remove contributions from 
both the primary CMB anisotropy and the kinematic SZ effect, as has been attempted in 
the analysis of MITO measurements of the Coma cluster (DePetris et al. 2002). 

The ratio of signals in two different photometric channels i and j is 



GiAQ\iIo (e*-ir{f dT 


'0fi(x) 


-P + R(x,9,pj\}ei( 


v)dv 


GiMlli SrT0r{fdT 


'oh(x) 


-P + R(x,e,pj\}ej 


'v)dv 



(4) 



The main dependence on the cluster properties in y cancels out when f3 is negligible. Multi- 
frequency observations that include measurements at the crossover frequency (e.g., MITO) 
afford effective separation of the thermal and kinematic components, exploiting their very 
different spectral shapes. This was demonstrated in the analysis of MITO measurements 
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of the Coma cluster (De Petris et al. 2002). Since the above ratio depends weakly on 
the cluster velocity, the residual uncertainty due to velocities of even ±500 km/s can be 
ignored in comparison with other errors. The ratio is also weakly dependent on the gas 
temperature at a level which we found to correspond to ~ 1% uncertainty in the estimation 
of the CMB temperature (for a typical observational error in T e ). Moreover, the uncertainty 
associated with the absolute calibration, Gi, is largely removed once we fit data from several 
photometric channels, as long as they are calibrated with a source with a known spectrum 
(e.g., a planet) even if its absolute calibration is uncertain. Only relative uncertainties among 
the various spectral channels are important; these include differences in angular response 
and in atmospheric transmittance. A standard blackbody source with a precisely calibrated 
temperature is therefore not required. We expect these considerations to imply that the 
precision of CMB temperature measurements via the SZ effect will not be appreciably affected 
by most of the known systematic errors. The level of precision in the measurement of T(z) 
is limited largely by other observational uncertainties, as discussed below. 

The responsivity Gi of each channel is usually determined from detailed observations 
of very well measured sources such as planets (mostly Jupiter or Saturn) and the Moon. 
While the temperature uncertainty of these sources can be as large as 10%, their spectra are 
relatively well known. For a source at a temperature of T s and with a throughput £1$, the 
signal in the Rayleigh- Jeans part of the spectrum is 



Since we are interested in the ratio of the signals in two channels, Ts drops out from the final 
expression, and the uncertainty in the ratio of values of the brightness temperature in the 
two channels depends only on the emissivity of the planet, rj(is), convolved with the spectral 
efficiencies of the channels, e«(z/). Note that the throughputs of individual channels, AQ\i, 
are usually made almost equal by an appropriate choice of the optical layout. 

The main uncertainty in the MITO measurements is due to imprecise knowledge of the 
bandwidth and transmittance of the filters. Therefore, we have precisely measured the total 
efficiency of our photometer by means of a lamellar grating interferometer. The response of 
the photometer is measured when it is illuminated with a laboratory blackbody at different 
temperatures ranging from liquid N 2 to room temperature. The effect of the uncertainty in 
our bandwidths has two main consequences. The first is the error in the ratio of responsivities 
of channel % and channel j, Gi/Gj. We estimate this error by computing the ratio of expected 
signals evaluated by convolving Moon and Jupiter spectra with the spectra of our filters, and 
taking into account spectroscopic uncertainties; we obtain a final error estimate which is less 
than 1%. The second consequence of bandwidth uncertainty enters in the integrations over 
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the SZ spectral functions in eq.[4]. However, this latter error is smaller than the former and 
is practically negligible. 

Also quite precisely measured were the angular responses of the telescope and the pho- 
tometer in the four channels. We have studied differences in the response of the four channels 
to both a point source, such as Jupiter, and an extended source, such as the Moon, as they 
cross the field of view. The optical layout of our photometer is such that all four detectors 
observe the same sky region. The differences among channels are less than 1% for extended 
sources. Finally, we have also studied the change in relative efficiency of our observations 
under different atmospheric conditions. The simulations show a change of ~ 1 mK in the 
estimated CMB temperature if the water vapor content changes from 0.5 mm to 1.5 mm, as 
determined by convolving theoretical atmospheric spectra with our filters. The uncertainties 
in determining the water vapor content, related to our procedure of measuring the atmo- 
spheric transmittance through a sky-dip technique, are at a level of 30%. However, when 
the effect of this error is propagated through our procedure and analysis we estimate that it 
amounts to a negligible error in the final result for the CMB temperature. 

In conclusion, systematic errors contribute less than 3% to the observed MITO signals, 
and are quite negligible with respect to the 10% contribution associated with the residual 
noise of the four detectors. For the OVRO measurement, we use the total error as specified 
by Mason et al. (2001); the statistical weight of this low frequency measurement is such that 
it only contributes ~ 12% to the final results. 

Our assessment of the errors in the SuZIE and OVRO & BIMA measurements of A2163 
(LaRoque et al. 2000) is far less certain than those in the MITO measurements, since we 
do not know the exact spectral responses of their filters and the atmospheric conditions 
during each of these observations. We have used Gaussian profiles for the SuZIE filters, with 
peak frequencies and FWHM values as given by Holzapfel et al. (1997a). For the OVRO & 
BIMA data point, we also take a Gaussian profile with a peak frequency at 30 GHz and a 1 
GHz bandwidth. We have employed data provided by SuZIE, for which a correction due to 
thermal dust emission has been taken into account. The Coma and A2163 data are collected 
in Table 1 and Table 2. Finally, the value we used for the electron temperature is (8.25±0.10) 
keV, as measured with XMM in the central region of Coma (Arnaud et al. 2001); the value 
for A2163 is 12.41? J keV (Holzapfel et al. 1997b). 

In both sets of cluster observations, we have minimized the difference between theoretical 
and experimental ratios (properly weighted by statistical errors) with T as a free parameter. 
The results are T Coma (0) = 2.7261q;qq4 K which, when multiplied by (1 + z) with z = 
0.0231 ± 0.0017, yields 

7b oma = 2.789l8;8™K. (6) 
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The main uncertainty is in the SZ observations, with only a small error (7 x 10 4 K) due 
to the uncertainty (0.1 keV) in the electron temperature. For A2163 we obtain 7^2163(0) = 
2.807l8$jf K , and when multiplied by (1 + z) with z = 0.203 ± 0.002, yields 

T A2 i63 = 3.377lg;X K. (7) 

The larger error in T e (~ 2 keV) translates to an uncertainty of 1.5 x 10~ 2 K in the deduced 
CMB temperature at the redshift of A2163. The mean of the above two values for T(0) is 
(T) = 2.7661^1 K > in g° od agreement with the COBE value. 

Using the COBE value of the temperature and the two values deduced here for T(z) at 
the redshifts of Coma and A2163, we fit the three data points by the relation T(z) = T(0)(1 + 
z) l ~ a \ the best fit is shown by the dashed line in Figure I. It corresponds to a = — O.I6I032 
at 95% confidence level (CL). With the alternative scaling T(z) = T(0)[1 + (1 + d)z], we 
obtain a best fit value of d = 0.17 ± 0.36 at the 95% CL. 

These first SZ results for T(z) are clearly consistent with the standard relation, although 
the most probable values of both a and d indicate a slightly stronger z dependence. Our 
values for these two parameters are quite similar to those deduced by LoSecco & Mathews 
(2001) from measurements of microwave transitions. Their results - based on two firm 
detections at redshifts 2.338 and 3.025 - are a = -0.05 ± 0.13 and d = 0.10 ± 0.28 (at 
95% CL). Thus, our SZ-based values already provide the same level of precision even though 
the two clusters are at much lower redshifts (z < 0.2). SZ measurements of higher redshift 
clusters will clearly provide a preferred alternative to the atomic and molecular lines method. 
If we fit the above relations to the combined data sets (consisting of all five temperatures), 
we obtain a = -0.09 ± 0.20 and d = 0.14 ± 0.28 (at 95% CL). Finally, Molaro et al. (2002) 
have recently questioned the validity of the temperature measurement at z — 2.338 (by 
LoSecco & Mathews 2001), and have revised the value at z = 3.025 to 12.1+3; 2 K. We have 
thus repeated the fit excluding the fourth data point and using the latter value for the fifth; 
the results from this fit are essentially unchanged: a = —0.11 ± 0.22, and d = 0.16 ± 0.32. 



4. Discussion 

The method employed here to measure T(z) can potentially yield very precise values 
which will tightly constrain alternative models for the functional scaling of the CMB temper- 
ature with redshift, and thereby provide a strong test of non-standard cosmo logical models. 
In order to improve the results reported here, multi-frequency measurements of the SZ effect 
in a significantly larger number (~ 20) of clusters are needed. The clusters should include 
nearby ones in order to better understand and control systematic errors. Exactly such ob- 
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servations are planned for the next 2-3 years with the upgraded MITO project (Lamagna et 
al. 2002), and the stratospheric BOOST experiment (Lubin P., private communication) and 
OLIMPO experiment (Masi et al. 2003). These experiments will employ sensitive bolometric 
arrays at four frequency bands with high spatial resolution. 

As detector noise is reduced, the uncertainty in the electron temperature may contribute 
the dominant relative error. For example, an uncertainty of ~ 2 keV - as in the case of 
Abell 2163 - becomes dominant when detector noise is reduced to a level ten times lower 
than the values used here. Therefore, it will be necessary to select clusters for which the 
gas temperature is precisely (to a level of ~ 0.1 keV) known. With the large number of 
clusters that have already been - or will be - sensitively observed by the XMM and Chandra 
satellites, it should be possible to optimally select the SZ cluster sample. 

With the currently achievable level of precision in intracluster gas temperature measure- 
ments, the present technique makes it possible to determine a and d with an uncertainty that 
is not less than 0.03 on both parameters, even with ideal (i.e., extremely high sensitivity) 
SZ measurements. The projected sensitivities of the future Planck and Herschel satellites 
will enable reduction of the overall error in the values of these two parameters by more than 
an order of magnitude. In addition to limits on alternative CMB temperature evolution 
models, such a level of precision will open new possibilities of testing the time variability of 
physical constants. For example, it will be feasible to measure a possible variation of the fine 
structure constant, thereby providing an alternative to the current principal method which 
is based on spectroscopic measurements of fine structure lines in quasar absorption spectra. 
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Fig. 1. — Black diamonds mark values of the CMB temperature obtained by COBE/FIRAS 
and SZ measurements. The two values reported by LoSecco et al. (2001) are shown as white 
squares, and the white diamond is the value reported by Molaro et al. (2002). The solid line 
is (l + z) scaling law as predicted in the standard model. The dashed line shows the best fit to 
the alternative scaling law with a determined by the COBE/FIRAS and SZ measurements, 
whereas dot-dashed lines are the ±lcr values of a. The inset shows a magnified view of the 
region z < 0.25. 



Table I. SZ signals from COMA used in the present analysis. 



Experiment Frequency [GHz] Bandwidth [GHz] AI [MJsr- 1 ] AT \jiK] 



OVRO 32.0 13.0 -0.0159 ±0.0028 -520 ± 93 

MITOl 143 30 -0.068 ±0.014 -179.3 ±37.8 

MIT02 214 30 -0.016 ± 0.039 -33.4 ±81.2 

MIT03 272 32 0.075 ± 0.015 169.8 ±35.1 
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Table 2. SZ measurements of A2163 used in the present analysis. 



Experiment 


Frequency [GHz] 


Bandwidth [GHz] 


AI [MJsr- 1 ] 


AT [fiK] 


OVRO+BIMA 


30 


1 


-0.048 ± 0.006 


-1777 ±222 


SuZIEl 


141.59 


12.67 


-0.380 ±0.037 


-1011.3 ±98.0 


SuZIE2 


216.71 


14.96 


-0.103 ±0.077 


-213.0 ± 159.3 


SuZIE3 


268.54 


25.66 


0.295 ±0.105 


662.2 ±235.7 



